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Abstract. Wc consider the Rcnyi entanglement entropies in the one- 
dimensional XX spin-chains with open boundary conditions in the presence of 
a magnetic field. In the case of a semi-infinite system and a block starting from 
the boundary, we derive rigorously the asymptotic behavior for large block sizes 
on the basis of a recent mathematical theorem for the determinant of Tocplitz 
plus Hankel matrices. We conjecture a generalized Fishcr-Hartwig form for the 
corrections to the asymptotic behavior of this determinant that allows the exact 
characterization of the corrections to the scaling at order o(^ _1 ) for any n. By 
combining these results with conformal field theory arguments, we derive exact 
expressions also in finite chains with open boundary conditions and in the case 
when the block is detached from the boundary. 
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1. Introduction 

The interest in quantifying the entanglement in the ground state of extended quantum 
systems has risen sharply in the last decade [1]. Remarkably, the results of many 
investigations allowed a deeper and more precise characterization of many-body 
systems. Furthermore, many surprising connections between fields and techniques 
apparently disconnected emerged. 

In this paper we will consider the Renyi entanglement entropies in the XX chain 
with open boundary conditions (OBC). As we will review, the leading asymptotic 
behavior of the entanglement entropy can be deduced directly from known results in 
conformal field theory (CFT) joined with available exact calculations for the chain with 
periodic boundary conditions (PBC). However, the open chain presents subleading 
oscillatory corrections to the scaling whose first observation dates back to 2006 [2] 
and that until now resisted to an analytic computation. This study provides a new 
and unexpected by-product that, for some aspects, is even more interesting that the 
main result itself. Indeed, in order to arrive to an analytic result for the entanglement 
entropy in systems with OBC, we faced the problem of calculating leading and 
subleading behavior of the determinants of matrices that in mathematical literature 
are known as Toeplitz plus Hankel (i.e. they are composed of a part that depends only 
on the difference between row and column indexes and another depending only on their 
sum). Toeplitz matrices have a very long history, culminating with the Fisher-Hartwig 
(FH) conjecture [3]. This conjecture has been proved (in some particular cases) only 
many years after its formulation by Basor [4]. The interest in the corrections to this 
formula leaded to a generalization known as generalized FH conjecture [5] that has 
not yet been proved. This formula has been fundamental to provide the corrections to 
the scaling for the entanglement entropy in systems with PBC [6, 7]. When moving 
from PBC to OBC, we move from Toeplitz matrices to Toeplitz plus Hankel ones, 
that is a brand new field of mathematics. The formula generalizing FH has been 
proved very recently [8] (see also [9]), and there is no conjecture for the subleading 
terms. Putting together the ingredients of the generalized FH and the recent results 
for Toeplitz plus Hankel, we conjecture a generalized FH formula for Toeplitz plus 
Hankel matrices, that we use to determine analytically the corrections to the scaling 
of the entanglement entropy for OBC. 

1.1. Entanglement entropy, conformal field theory, and boundaries 

Let us consider an infinite one-dimensional critical system whose scaling limit is 
described by a CFT of central charge c, and a partition into a finite block of 
length I and the remainder. The entanglement entropy (the von Neumann entropy 
— Tr pa log pa of the reduced density matrix pa) for t much larger than the short- 
distance cutoff a is asymptotically [10, 11, 12, 13] 



where c'j is a non- universal additive constant. A more general measure of the quantum 
entanglement is provided by the Renyi entropies, that in the conformal situation of 
before scale like [12, 13] 



Si = - 



Tr p A In p A ^ x In - + c[ 



a 



(1) 
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Figure 1. The two block configurations we consider. Top: a block of length £ 
starting from the boundary. Bottom: a block of length i at distance l§ from the 
boundary. 



where c' n are other non-universal constants. For n = 1, Eq. (2) reproduces Si, but 
the knowledge of S n for different n contains more information and characterizes the 
full spectrum of non-zero eigenvalues of pa [14], that is fundamental to understand 
the scaling of some numerical algorithms based on matrix product states [15]. 

For semi-infinite systems, Renyi entropies depend on where the block is placed 
because of the breaking of translational invariance. When the block starts from the 
beginning of the semi- infinite chain (see Fig. 1), very general results are available from 
boundary CFT. The Renyi entropies behave like [12, 13] 

s»^fi + rW + ^> ( 3 ) 



12 

The constants c' n are non-universal, but their value is related to c' n for systems with 
PBC by the universal relation [12, 16, 2] 

4-y =logff, (4) 

where logg is the boundary entropy, first discussed by Affleck and Ludwig [17]. g 
depends only on the boundary CFT and its value is known in the simplest cases. For 
the XX model with OBC considered here, we have g — 1 [17]. 

When the block of length £ placed at a distance Iq from the boundary (see Fig. 
1), the situation is more complicated. In this case, global conformal invariance only 
gives the general scaling of S n as 

Sn = 12 ^ + n ) bg Wmi + to) + 2C » + -T^n~ ' (5) 
where x is the four-point ratio 

e 



" (24 + 1) 2 ■ (6) 
This formula can be readily obtained from the entanglement of two disjoint intervals 
in an infinite system [18, 19], where the points corresponding to the second interval are 
the mirror images (with respect to the boundary) of the actual interval. The function 
F n {x) (normalized such that ^(0) = 1) depends on the full operator content of the 
theory and must be calculated case by case. The entanglement of two blocks attracted 
in the last year an enormous interest [18, 19, 20, 21, 22, 23] because of its ability to 
detect the full operator content of a CFT. However, all the results presented so far are 
for periodic systems and still little attention has been devoted to the boundary case 
that we address here in a specific case (that we will see to turn out to be trivial, after 
the corrections to the scaling have been properly taken into account). 
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1.2. Beyond CFT 

When S n (t) is computed numerically, it has been observed that the asymptotic CFT 
result is obscured by large, and often oscillatory, corrections to the scaling [2, 6, 24]. In 
Rcf. [6], on the basis of both exact and numerical results, it has been argued that these 
corrections arc in fact universal and encode information about the underlying CFT 
beyond what is captured by the central charge alone. More precisely, they give access 
to the scaling dimensions of some of the most relevant operators of the underlying 
CFT. This conjecture of Rcf. [6] has been recently confirmed by using perturbed CFT 
arguments [25] (and generalized to gapped chains as well [26]). For a Luttingcr liquid, 
the proposed scaling form of S n is [6, 25] 

S n = S° FT + f n cos{2k F l)r 2K ' n , (7) 

where K is the scaling dimension of a relevant operator (in general the oscillating factor 
can be different from cos(2fcp£) or even be absent, as it happens for the Ising model 
[7, 27]). The constant /„ is a non-universal quantity that has been determined exactly 
only for XX and Ising models [6, 7]. It is worth mentioning that in the known PBC 
cases /i turns out to be zero, i.e. for the entanglement entropy there are no unusual 
corrections. This is also the case for all the numerical computations presented so far, 
but this fact still lacks of a general proof. 

On the basis of universality [6] and directly by CFT [25] , it has been argued that 
in the case of OBC, the exponent governing the corrections is half of the PBC one 
(i.e. K/n replaces 2K/n in Eq. (7)), as compatible with all numerical computations 
available [2, 28, 29]. It is important to mention that with OBC, unusual corrections are 
also present for n — >• 1 [2]. In this manuscript, we present a first analytic computation 
of these corrections to the scaling for systems with OBC. 

We anticipate here the main result of the manuscript for the Renyi entanglement 
entropy of a block starting from the boundary of a semi-infinite system 

S n (l) = i (l + ^ In [2 {21 + 1) | sinM] + ^ 

+ 2Sln[ ^ +1) W + Dl sink F \]-^ T -$±4\ + o(rV») , (8) 

that is compatible with the general Luttinger liquid prediction [6, 25] with K = 1. 
Notice that away from half-filling (hp = n/2), the oscillations have different forms 
compared to the PBC case. While this formula is correct only to order o(£ _li/ ™), in 
the following we will present the full expansion up to the order o(^ _1 ), for any finite 
value of n. In the case of n —¥ oo, the corrections become logarithmic and are also 
exactly calculated. CFT is also used to infer exact analytic formulas for finite systems 
with two open boundaries. 

1.3. Organization of the manuscript 

The remainder of this paper is organized as follows. We first present known results for 
the XX model and the calculation of corrections to the scaling with PBC in Sec. 2. 
In Sec. 3 we derive the asymptotic behavior of Renyi entropies, while in the Sec. 4 we 
calculate analytically the first correction to the scaling and a family of the sublcading 
ones that, independently of n, give the correct behavior at o(^ _1 ). In Sec. 5 we 
generalize the results to finite systems. In Sec. 6 we consider the case of a block 
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disconnected from the boundary. Finally in Sec. 7, we summarize our main results 
and discuss problems deserving further investigations. 

2. Entanglement entropy in the XX model 

The Hamiltonian of the XX model for a semi-infinite chain with an open boundary is 

1=1 

where a^' y ' z are the Pauli matrices at site I. The Jordan- Wigner transformation 



\m<l J 

maps this model to a quadratic Hamiltonian of spinlcss fcrmions 

H = - E c ? c '+i + c li c; + 2/1 - • ( n ) 

Here h represents the chemical potential for the spinless fermions c/, which satisfy 
canonical anti-commutation relations {q,cJ„} = 5^ m . The Hamiltonian (11) is 
diagonal in momentum space and for \h\ < 1 the ground-state is a partially filled 
Fermi sea with Fermi-momentum 

k F = arccos \h\. (12) 

In the following we will always assume that \h\ < 1 so that we are dealing with a 
gapless theory. 

Using Wick theorem, the reduced density matrix of a block A = [£ + 1,£ + £] 
composed of £ contiguous sites in the ground state of the Hamiltonian (9) can be 
written as 

p A = detC exp ( [HC- 1 - 1)] jt c] Cl I , (13) 

where the correlation matrix C has matrix elements defined by 

C nm = (cLc„) (14) 

In the ground-state of the open semi-infinite Hamiltonian (9), the elements of the 
correlation matrix are [31] 

_ sin(fc F (n - m)) sin(k F (n + m)) 

<^nm — / \ 7 . \ j (15) 

ir(n — m) Tr(n + m) 

and the PBC result is recovered when n, m — > oo while keeping their distance n — m, 
finite, i.e. far from the boundary as the physical intuition suggests. The first half 
of this expression (Toeplitz part) is the same as in systems with periodic boundary 
conditions, while the second part (Hankel type) is direct consequence of the non- 
translational invariant terms introduced by the boundary. 

As a real symmetric matrix, C can be diagonalized by an orthogonal 
transformation 

RCR T = S lm (l + u m )/2, (16) 
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and the eigenvalues depend both on £ and £q. The reduced density matrix pa is 
uncorrelated in the transformed basis, so that the Renyi entropies can be expressed 
in terms of the eigenvalues vi as 

£ 1 

S n {£ J) = y2e n (v l ), with e n (x) = - In 

z — ' 1 — n 

1=1 

More details about this procedure can be found in, e.g., Rcfs. [11, 30, 31]. The 
above construction refers to the block entanglement of fermionic degrees of freedom. 
However, the Jordan- Wigner transformation, although non local, mixes only spins 
inside the block. As well known, this ceases to be the case when two or more disjoint 
intervals are considered [21, 32] and other techniques need to be employed [33] in 
order to recover CFT predictions [18, 19, 21]. Furthermore, also in the case of XX 
chains with different boundary conditions (e.g. fixed) the Jordan- Wigner string would 
spoil the correspondence between spins and fermions for an interval detached from the 
boundary (i.e. £q ^ 0). It is a peculiarity of OBC that the reduced density matrix of 
any interval at any distance from the boundary is the same for spins and fermions. 

The representation (17) is particularly convenient for numerical computations: 
the eigenvalues v m of the I x £ correlation matrix C are determined by standard 
linear algebra methods and S n (£) is then computed using Eq. (17). All the numerical 
computations presented in the following have been obtained in this way. 

The sum in Eq. (17) can be put in the form of an integral on the complex plane 
[34], introducing the determinant 

D/(A) = det ((A + 1)1 - 2C) = det(G) . (18) 

In the eigenbasis of C the determinant is simply a polynomial of degree £ in A with 
zeros = 1, . . . ,£}, i.e. 

i 

D t {X) = H{X-uj). (19) 

This implies that the Renyi entropies have the integral representation 

S n (£) = ±f dX e n (X) dh ^ , (20) 

where the contour of integration encircles the segment [—1,1]. In the PBC case 
and in the thermodynamic limit (L — > oo), Fishcr-Hartwig conjecture allows to 
obtain the asymptotic large £ behavior of S n (£) [34] . The generalized Fisher-Hartwig 
conjecture permits the computation of all harmonic corrections [6, 7], while non- 
harmonic corrections can be computed only exploiting random matrices techniques 
[7]. In next subsection, we report the Fisher-Hartwig approach to PBC, in order to 
fix the notation and to understand the needed ingredients for OBC. 




2.1. The asymptotic result for periodic boundary condition. 
For PBC, the correlation matrix in the limit L — > oo is 
_ sm(k F (n-m)) 

^nm — 7 n • \^1) 

TT(n — to) 

The matrix G in Eq. (18) is a £ x I Toeplitz matrix, i.e. its elements depend only on the 
difference between row and column indices Gjk = g^-k- I n this case, it is possible to 
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use the (generalized) Fisher-Hartwig conjecture to calculate the asymptotic behavior 
of Di(X) and hence the Renyi entropies [34, 7]. The standard FH calculation proceeds 
as follows. We define the symbol of the Toeplitz matrix the Fourier transform g(9) of 

m 



9l 



dO 
2^ 

that in our case takes the form 
'A + l 



e m 9(0), 



(22) 



9(0) 



(23) 



9 £ [k F ,2-n- - k F ] 
v A- 1 9 £ [0,k F ] U [2tt - k Fl 2n] . 

On the interval [0, 27r] the function g(6) has two discontinuities at 9i = k F and 
9 2 = 2tt — k F . In order to employ the Fisher-Hartwig conjecture one needs to express 
g(6) in the form 

R 

9(d) = fifi) ]J e^ 9 - 8 "-^ 9 ^)] ( 2 _ 2 cos(9 - 9 r )) a " , (24) 



where R is an integer, a r , b r and 9 r are constants and f(9) is a smooth function with 
winding number zero. The Fisher-Hartwig conjecture then states that the large-^ 
asymptotic behavior of the Toeplitz determinant is 



D t - F[f(9)} 




(25) 



where F[f(9)] — exp(^ J Q 27r d9\nf(9)) and E is a known function of f(6), a r , b r , and 
9 r . In our case it is straightforward to express the symbol in the canonical form (24). 
As g(9) has two discontinuities in [0, 27r) we have R — 2. By comparison, we have 

ai,2 =0, 

b 2 = -b 1 =f}\ + m, 

f(0) = fo = (A + l) e - 2ib2kF = (A + l) er ^k F m e ^k F p x ^ ( 2 g) 

where m is an arbitrary integer number, that labels the different inequivalent 
representations of the symbol g(0), see [4], and 



2m 



A + l 



A 



A + l 

— 7r < arg < 7T . 

A — 1 



(27) 



Jin and Korepin employed the Fisher-Hartwig conjecture for the m = representation 
and obtained the following result for the large-£ asymptotics of Df (\) [34] 



(A + l) 



A + l 
A- 1 



(2*| sinM)- 2 ^ A )G 2 (l + /3 A )G 2 (1 - f3 x ) , (28) 



where G(x) is the Barnes G- function [36]. Inserting (28) into (20) and carrying out 
the integral leads to the result for the asymptotic behaviour of the Renyi entropy 



1 + 



1 



1 



n J 



ln(2*|sinfc F |) +E n . 



where the constant E n has the integral representation 



Ey, 



1 



00 dt 

T 



l 



l 



l 



nsmht/n sinhty sinht 6 



1 



(29) 



(30) 
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However, when the symbol has several inequivalent representations as in our case, 
the generalized Fisher-Hartwig conjecture (gFHC) [4] applies and one has to sum over 
all these representations as 

Di(X) ~ J2(f (m)) e r^r=^ m » 2 E(m) , (31) 

m 

where all the various FH constants /o, h, and E depend on m, as shown in Eq. (26). 
Then the full result of the generalized Fisher-Hartwig conjecture for the Tocplitz 
determinant takes the form [7] 

D^iX + iyf^) * ]T(2^| sin fc F |)-2(™+^)V2*M 

x [G(m + 1 + X )G(1 - ro - /3 A )] 2 , (32) 
leading after long calculations for the integral (20) to [7] 

1 n l L \2 2TJ- 

We mention that the gFH conjecture has also other applications in physics as for 
example those reported in Refs. [37, 38]. 



+ o(r 2/n ). (33) 



3. Fisher-Hartwig like conjecture for the semi-infinite chain 

We want to use a generalization of the FH conjecture to obtain the determinant in Eq. 
(18) for large I. The matrix G is not of the Toeplitz form because of the second term 
in Cnm in Eq. (15). Let us consider a block at distance £q from the boundary (i.e. 
starting from the site £o + 1). To fix the notation, the matrix G is an (l+l) x (£ + 1) 
matrix that is a sum of Toeplitz and Hankel matrices with elements of the form 

G nm = 9n-m ~ 9n+m+2£ + 2 , 71, 771 = 0, 1, . . . i - 1 . (34) 

Here gi is the same as in Eq. (22). We are not aware of any possible generalization 
(conjectured or proved) of the FH formula to the case of arbitrary £ . However, very 
recently a theorem has been proved [8] for some Toeplitz+Hankel matrices, among 
which those of the form 

a i-j ~ a i+j+2 j (35) 

that applies to our case with £q = 0, i.e. when the block starts from the boundary. 

For the case of a piecewise constant symbol as the one of our interest (cf. Eq. 
(23)), this theorem was known from longer time [9]. The result in Refs. [8, 9] has 
a structure similar to the FH formula, with more constants aj , bj corresponding to 
boundary terms. Such a general formula is not very illuminating and too long to be 
written in its full glory. We remand the interested reader to the original reference [8] . 
We specialize this formula to the symbol in Eq. (23), and we obtain the asymptotic 
behavior of the determinant Di(X) as 

D e (X) ~ e ^- k ^ [(A + 1) __f ] \u\ sm{k F )\)-' f 'G(l - /3)G(1 + 0) , (36) 

where (3 has been defined in Eq. (27). The result is very similar to the PBC case. The 
main differences are in the halving of the exponent of £, in the factor 4 instead of 2 
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multiplying I (both with a clear physical interpretation) , in the absence of the square 
in the Barnes G function, and in the phase shift in front. 

We can now proceed to calculate the integral (20) giving the asymptotic Renyi 
entropies. First we notice that the phase shift e l ^~ kF ^ in Eq. (41) gives a vanishing 
contribution to the entropies, indeed 

V2 7T / 2-Ki J 1 - A 2 y J 

The remaining part of the integral parallels that for PBC, giving 

S n (£) = ^{ 1 + l) sinfc F |) + ^ , (38) 

where E n is defined in Eq. (30). 

This result for S n is exactly what expected from CFT, as reviewed in the 
introduction. Both the leading logarithmic term and the subleading constant one are 
in agreement with Eq. (3) with ln g = 0. There is then no new physical information 
in this expression. However, the present result is based on a mathematical theorem 
and so it provides a rigorous confirmation of a general CFT result in a specific lattice 
model. The new physical results of this paper are in the following sections. 



4. Generalized conjecture for the oscillating corrections to the scaling 

We consider now the corrections to the scaling to the asymptotic result derived in 
the previous section. Analogously to the PBC case, we need a generalized Fisher- 
Hartwig formula that applies to the case of Toeplitz+Hankel matrices. To the best of 
our knowledge, this formula does not exist neither conjectured or proved. However, 
armed with the general ideas presented above is not difficult to give our own conjecture 
for such a determinant: we should only sum on the several different inequivalent 
representations of the symbol. 

For the symbol in Eq. (23), we conjecture the following formula 



£>/(A)~(A + l)< 



A + l 

A- 1 



e if(/3+"0 



sinfct 



x e -2*fcH£+m)(£+i/2) G(m + i + j 8 A )G(i - m - /3 A ). (39) 

Most of the above formula is inspired to the gFH Eq. (32) and adapted to the present 
case. However, the factor 1/2 as an additive constant to I has been introduced without 
any mathematical reason. This factor 1/2 gives an analytic (i.e. non-harmonic) 
correction to Dg(\) and we introduced it to reproduce accurately the numerical data. 
In Fig. 2 we report the ratio of the numerically calculated Di(X) asymptotic value 
with and without the additional 1/2, showing that the former converges faster. The 
data for the resulting entanglement entropy S2 (t) (reported in the right panel of Fig. 
2) show even more clearly the importance of this factor. In any case, we conjectured 
this gFH formula and in doing so we prefer to conjecture an expression that reproduces 
numerical data as accurately as possible. A full justification of the factor 1/2 could 
be obtained by generalizing the random matrix results for PBC in Rcf. [35], but this 
is beyond our knowledge. 
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Figure 2. Left: The ratio R = D t (X) / Df' (\) for A = 1 + i and k F = ^JV[2. 
Right: The entanglement entropy for kp = y/n/2. Both graphs show the 
importance of the additive factor 1/2. 



The leading corrections to the scaling is obtained by summing only the three 
modes m = — 1, 0, 1, obtaining 



Dp ~ \ l + ie 



-ikf 



-2ik F l 



L 



^D?\l + * e (X)), 
where we isolated the leading term 



-1-2/3 r(i + /?) _ iclkFr2lkFlT7l+20 T^~fi) \ 



T<J3) I 



D (0) = e i(f-fc,)/3 



(A + 1) lX^l) " 1 ^ 2 G(l-^)G(l + /3) 
we defined ^(A), and for OBC is defined as 
L k = 2(2£+l)\sm(k F )\. 

We define 



d n (£) = S n (£) - S° n (£) , 



with 



c(o) 



12 V n 



In [2(2^ + 1) |sinfc F |] + 



-En 



(40) 

(41) 

(42) 
(43) 

(44) 



that includes the additive 1/2 factor compared to the leading term in Eq. (38). For 
large L k we have 

1 

2ni T dX 



d n (£) 



dX e n (X) = M? dX en{X) ^XT 



(45) 



The contour integral can be written as the sum of two contributions infinitesimally 
above and below the interval [— 1, 1] respectively, i.e. 



d n (£) 



1 


■ pl+ie 


pl—ie " 


2~7Ti 


J-l+ie 


l-l-ie. 



dX e„(A) 



d*i(X) 
dX ' 



(46) 
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This shows that we only require the discontinuity across the branch cut. The only 
discontinuous function is j3\, which for —1 < x < 1 behaves as 



Px±te = -iw(x) T 2 ' 
We now change variables from A to w 



with 



( \ 1 i 1 

w(x) = — In 

w 2tt 1 



A = tanh(7rw) 



-co < w < oo. 



(47) 



(48) 



We have 

-1-20 r(i 



■P) 



-1+2/ ^(1-/3) 



[3 — — iw— \ 



„, r 2iw 
— L k 



-iw-\-4 



where we have dropped terms of order 0(L k 2 ) compared to the leading ones and we 
have defined 



iw) 



va + iw) 



Integrating by parts and using 
d 



, e„(tanh(7rw)) = - 
dw 1 



irn 



(tanh(n7rw) — tanh(7rw)) . 



(49) 



(50) 



we arrive at 



in f 00 

d n (£) ~ — -. r / dw(tanh(7rw) — tanh(n7T'u;)) 

2(1 - n) J_ QO 



ie- lkF e- 2lkFl Ll iw 



+ ie ikF e 2ik F e L -: 



,r(|+^)i 



(51) 



r(| + iw) 1 " " fe v{\ -iw) 

For large t the leading contribution to the integral arises from the poles closest to the 
real axis. These are located at wq — i/2n (wq = —i/2n) for the first (second) term. 
Evaluating their contributions to the integral gives 

-i/nT(- + —) 



d n (£) 



2sin[fc F (2^+l) 



2(2£+ l)|sinfc F | 



1 V 2 2n I 



(52) 



In Fig. 3 we report the numerical calculated S n (l) for a semi-infinite chain and 
for small values of n = 1,2,3. For these values of n, the inclusion of only the first 
correction to the scaling (as in Eq. (8)) is enough to describe very accurately S n {t) 
even for relatively small values of £. The figure shows the correctness also of the fc F 
dependence of the correction, that is the most important difference compared to PBC. 
Notice that without the factor 1/2 introduced by hand in Eq. (39), the correction (52) 
would be inadequate at order l/£, as anticipated in Fig. 2. 



Subleading corrections 

Eq. (33) describes the asymptotic behavior in the limit Lk — > oo with n fixed. 
It provides a good approximation for large, finite £ as long as ln(Lfc) » n. For 
practical purposes it is useful to know the corrections to S n (£) for large I but ln(L^) 
not necessarily much larger than n. In this regime there are two main sources of 
corrections to (33). The integral (50) is no longer dominated by the poles closest 
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Figure 3. S n (£) in a semi-infinite XX chain for n = 1, 2, 3 and different values of 
kp specified in the captions. The exact numerical results are compared with the 
asymptotic formula including only the first correction to the scaling. For these 
small values of n, the agreement is excellent. Notice the incommensurability 
effects when kp is not a fraction of it (lower graphs). 



to the real axis and contributions from further poles need to be included. These 
give rise to corrections proportional to L k q ^ n , with q integer. Furthermore, terms 
in the expansion of the logarithm in Eq. (45) need to be taken into account. The 

corresponding contributions are proportional to e ±l2pfcF ^ with p = 2, 3, 

We now take both types of corrections into account. The following derivation is 
very similar to the one for PBC in Ref. [7]. We first consider the series expansion of 
the logarithm in Eq. (45). 

h[1+ *, w]= f; (-"' +1 (W . (53) 

P =i 

Recalling the explicit expression (40) for ^i(A) leads to a binomial sum 
(* £ (A)) P = (ze^e- 2 ^ e L- (1+2M c 0x - ie~ lk * e 2lk ^ L,- (1 - 2M c_p x ) P 

= J2 ^\p(-l)« e l ^( 2 ^P)e 2lfcF£ ( 2 «-rfL- p L^ 2( ^ 29)/3A C P7cl /3A , (54) 

where we have introduced the shorthand notation c« = T(l + /3)/T(—f3). 
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When calculating the discontinuity across the branch cut running from A = — 1 
to A = 1 all terms other than q = and q = p give rise to terms that are subleading 
in Lk- Hence we may approximate 

(\l>£(tanh(7rw) + ie)) p - (\Ev(tanh(7ra;) - ie)) p w 



j P e -2ik F (l+l/2)p L 2iwp ( p 

Plugging this into Eq. (45) 
~ (_i)p+i 



-1/2 



i y e 2ik F (l+l/2)p L -2iwp c p 



+iw-l/2 



(55) 



P 



2(1 



/C 



dw(tanh(7rw) — tanh(n7rw)) 



P =i 

x i p \e- 2l P k ^ e+1 /VLl lwp j p (w) + (-l)P+ 1 e 2ipk ^ e+1 ^Lf iwp Y(-w)] . (56) 

The integral is carried out by contour integration, taking the two terms in square 
brackets into account separately. The first (second) contribution has simple poles in 



the upper (lower) half plane at w q 
integer such that 2q — 1 ^ n, 3n, 5n, 



1 K 



where q is a positive 



" 2n \ w 1 " 2n 

for n 7^ 1. Contour integration then gives 



d n (i) 



(-1)p +1 (p{-K-2k F ) 



COS 



l - n ^ p V 2 

p,q=l 

where we defined the constants Q n , q as 

Qn.q 



2pk F e)L k Pi2q " 1 \Q n , q ) p + ..., (57) 



'2 



2a; 



•■2 2n ) 

In the sum over g the special values 2q — 1 7^ n, 3n, 5n 



(58) 



. are to be omitted for n 7^ 1. 
Eq. (57) shows that there are contributions to the Renyi entropies with oscillation 
frequencies that are arbitrary multiples of 2k p. 

At half-filling (kp = n/2) certain simplifications occur. For even £ we find 



d n (£) 



l-n 



(4*)— q„,i - {uy 



Qn,l 



Qn,l 



Qn,3 



while for odd £ we obtain 



d n (£) 



1 



(M)-^Q nA + (uy 



Ql 



+ m- 



Ql 



+ Q 



n,3 



+ 



(59) 



(60) 



that are of the same form as for PBC [7] with exponents that are halved. However, 
we stress that this is true only at half filling, while for generic fcp, Eq. (56) shows 
oscillations different from its PBC counterpart. 

In all the above analysis we have ignored contributions to the generalized Fisher- 
Hartwig conjecture with \m\ > 1. While these lead to oscillatory contributions with 
frequencies that are integer multiples of 2kp they arc suppressed by additional powers 
of and hence are subeading, even in the case where n is not small. 

In Fig. 4 we show the corrections d n {£) for n = 20 and 50 with kp — tt/6 and 
7r/2 respectively and their comparison with the asymptotic result Eq. (57). Step by 
step we take into account further terms in the asymptotic expression until we obtain 
a satisfying agreement with the numerical data. For n — 20, 13 terms in Eq. (57) arc 
enough to reproduce the data, while for n — 50 we need 41 terms to have the same 
accuracy. These numbers are larger than the corresponding ones for PBC [7] because 
the corrections in the present case have smaller exponents. 




Figure 4. Difference of the Renyi entanglement entropy with its asymptotic value 
for n = 20 (left) and 50 (right) with kp = tt/6 and tt/2. The exact numerical 
data are compared with Eq. (57) with increasing number of terms. For n = 50 we 
report curves with 1,3,5,9,15,25,41 terms, while for n = 20 we consider 1,2,3,5,8,13 
terms. Increasing the number of terms considered, the expansion (57) becomes 
accurate even for values of i as small as 2. 



4-2. The limit of large n 

It is apparent from (57) that the limit n — > oo deserves special attention. Soo(£) is 
known in the literature as single copy entanglement [39]. Here it is necessary to sum 
up an infinite number of contributions in order to extract the large- asymptotics. It 
also provides information on the behavior of S n (£) in the regime n 3> lnL^, 1. 
The following derivation parallels the one for PBC in Ref. [7]. 

In order to investigate the limit n — > oo we consider Eq. (56), but now first take 
the parameter n to infinity and then carry out the resulting integrals. This gives 

doo{£) ~ t. ^2 / dw(sgn(w) - tanh(7rw)) 

1 P =i p J -°° 

x i p J e -2^(f+i/2)p L 2™ P [7(w)]P + (_ 1 )P e 2ifc,(m/2)„ i -ait U p [ 7 (_ W )]P 

-2 lfcF (£+l/2) Plm / _ tanh (™)] L 2™P [7(w)]P 



(_l)f e 2ife^(^+V2)P Im / dw[l - tanh{nw)}L~ 2lwp [ 1 {-w)] p 



.(61) 



Using that the first singularity in the upper (lower) half plane occurs at w = i/2 
(w = —i/2) we deform the contours to run parallel to the real axis with imaginary 
parts z/4 and — i/4 respectively, i.e. for the first term we use 

1 ri/^ /■oo+i/4 

duu f(w) = / dw f(w) + / dw f(w). 

Jo Ji/4 

It is straightforward to show that the second integral contributes only to order 0(1/ Lf.) 
and does not give rise to logarithmic corrections. Hence the leading contribution is of 



Parity effects in the entanglement of XX chains with open boundary conditions 15 




k F =7L^ 



Figure 5. l/\doo(i)\ with d x (£) = S x (e) - S { £> {£) vs I in log-linear scale for 
kp = 7r/2. Corrections to the scaling are logarithmic and perfectly described by 
Eq. (65) that is represented as two straight lines for I even and odd respectively. 



the form 



dooO ~ E 
p=l 



{-if 



1/4 



e - 2l fc F (£+l/2)p Re / ^ _ . tan ( 7rz )) jL -2^ ( 7 ( iz ))P 



1/4 



(-l)f rp{-K-2k F ) 



p=i ^ 



1/4 



-2pfcf*) / dze- 2 ^ lnLfc ( 7 (iz)) p . (62) 



For large Lfc the dominant contribution to the integral is obtained by expanding 
(7(iz)) p in a power series around z = 

rl/4 



4M~2V^ L cos 



{-If fp(n-2k F ) 



p=l 

OO 



P 



2pfc^ / dze~^ plnijfc (l + 2pz*(l/2) + ...) 



2 2J cos ( — o — — - 2 P fc i^ 



0=1 



p 



1 g(l/2) 
2 In L fc 2 In 2 L fc 



(63) 



where 4'(^) is the digamma function ^(z) = r'(z)/r(z). 

As in the PBC case, at half-filling (k F = it/2) this result takes a particularly 
simple form 

1 IT 2 

(-l) p(f+1) _ ) lnL fc 6 

1 it 2 



1 



InL, 



E 

P =i 



£ odd 



(64) 



\nL k 12 



even . 



Summing some of the subleading terms in (63) to all orders in (ln(Lj,)) 1 leads to an 
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expression of the form 

7T 2 f 2 I odd , 

d °°* ~ 191 JhT A 1 V ^ 

l2ln(bL k ) 1—1 £even, 

where b = exp(— \P(l/2)) w 7.12429. To check this result, in Fig. 5, we report the 
exact numerical data (only for fop = ir/2) for l/|cZ 00 (£)| in log-linear scale, showing 
explicitly the logarithmic form of the corrections, described very precisely by Eq. (65). 



5. Finite systems 

The Hamiltonian of a finite XX chain with two open boundaries is 

i=i i=i 

As before, the Hamiltonian is diagonalized by a Jordan- Wigner transformation and 
Fourier transform. However, we are dealing with a finite chain at fixed magnetic field, 
that in the language of fermions is fixed chemical potential. In this case, the number of 
fcrmions is a non-continuous function of L, because it can assume only integer values. 
The number of fermions is Np = [(L + l) \ arccos h\\. There is an ambiguity in defining 
the Fermi momentum and we choose the definition 

k F = = y[(L + \ arccos h\\ , (67) 

i.e. k F is an integer multiple of n/L. This definition has the advantage to give 
kp = ir/2 at half-filling for L even. However, the correlation matrix C is better 
defined in terms of 

fc ^^LTT + 2(ZTi)- ™ 

Notice the following important properties: 1) In the limit L — > oo, kp = k' F : 2) If 
kp = 7r/2, then k' F = ir/2; 3) For L odd, the ground state is doubly degenerate 
and both kp and k' F cannot be equal to ir/2. With this definition, we can write the 
correlation matrix as 



C ! 



2(L + 1) 



sin(k' F (n — mj) sin(fc^(n + m)) 



sin (2{LTi)( n - m )) sin (2(fTi)( n + m )) 



(69) 



From this correlation matrix, it is straightforward to obtain numerical results for 
S n {t) also in finite systems. This analysis has been already done with considerable 
numerical accuracy in Refs. [2, 29] for the Von Neumann entanglement entropy and 
in Rcf. [28] for general n. However, the accurate results for the amplitudes of the 
corrections to the scaling were not compared with the theoretical predictions not 
available at that time. 

The modification of the leading term in Eq. (8) for finite systems is provided by 
conformal field theory [12]: the Renyi entropies are given by Eq. (8) where the length 
of the subsystem I is substituted by the chord distance L/irsm(ir£/L). However, we 
would like an expression that takes into account corrections to the scaling and that is 
accurate at order l/£, while we keep fixed the ratio IjL. This is beyond the predictive 
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Figure 6. Embedding of the open finite chain in a periodic one. The blue points 
correspond to the original chain, the red ones are the mirror symmetric sites, 
while the two yellow points at x = and x = L + 1 = — L — 1 arc "auxiliary 
sites" . The picture shows a block with £ = 2 in the open chain corresponding to 
2£ + 1 = 5 in the periodic one. 



power of CFT, but an intuitive argument to find a proper modification of Eq. (8) 
valid in finite size can be given, leading to the expression: 



S n (£) = 



1 

12 



1 



1 



In 



4(i + 1) ^ tt(21+1) 
2(L + l) 



■ sin 



sin kw 



En 

~2~ 



2sin[fc^(2^ + 1)] r4(L + 1) . tt(2£+1) 
1-n I n Sm 2(L + 1) 



sin fcp 



-Vnr(i + 



2n i 



r( 



2n' 



(70) 



The argument proceeds as follows. While we take the continuum limit from the spin- 
chain to the CFT, there is a well-known arbitrariness on the exact correspondence 
between the lattice sites and the coordinate on the continuum space. While for PBC, 
translational invariance guarantees that we can start the lattice in an arbitrary point, 
this is no longer true in the presence of boundaries. For a semi-infinite system, the 
exact result (8) suggests that the first site of the chain should be placed at position 
x = 1 in the continuum theory. Indeed, when building the mirror image (as usually 
done in boundary CFT), we have a mirror chain starting from —1 going up to — oo. 
This implies that an "auxiliary site" should be introduced at x — 0. In this way, we 
have an infinite chain with a block of length 2£+l, exactly as Eq. (8) suggests. When 
we move to a finite chain of length L, the mirror construction is graphically depicted 
in Fig. 6. We clearly have to add another auxiliary site at the other boundary to 
embed the open chain in a periodic one. The resulting length of the periodic chain 
is 2(L + 1). Thus this argument suggests that from the semi-infinite formula (8), 
we can obtain a finite-size ansatz by replacing 2£ + 1 with the modified chord length 
2< " L ^^ sin ^2(l+i) • ^ n doing so, we should also keep in mind that the prefactor of the 
correction s\n[k' F {2£ + 1)] is not scaling (as for PBC [6]) and there i should be left 
unchanged. All these ingredients lead to Eq. (70). 
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L=152 £ k'p=Ji/2 L/(L+1) L=153 e k' F =Jt/2 ML+1) 



Figure 7. (£) in a finite XX chain of length L for n = 1,2,3 and different values 
of kp. The exact numerical results for any i are compared with the asymptotic 
formula including the first correction to the scaling. The excellent agreement 
confirms the effectiveness of the finite-size scaling ansatz Eq. (70). 

In Fig. 7, we report numerical calculated S n (£) for finite systems of different 
lengths and for different values of kp- In all cases, for n = 1,2,3, when the first 
correction describes accurately the numerics for semi-infinite systems, we found perfect 
agreement between analytical and numerical results, confirming the validity of the non- 
rigorous argument reported above. Notice that when k' F is not a simple number, as in 
the two graphics in the bottom of Fig. 7, the numerical data (points) show apparently 
strange periodicity. When the asymptotic exact forms are plotted (continuous lines), 
it is clear that the periodicity is the correct one and the previous effect is only due to 
the the value of kp. 

Notice that at n = 1, the unusual correction in l~ x l n and the analytic one coming 
from expanding 21 + 1 in the leading term are of the same order 1 /£ and they have 
been disentangled in the past [2, 29] only because one is oscillating and the other 
is not. The presence of the non-oscillating 1/L term has been firstly observed and 
its analytic value guessed in Ref. [2] for n — 1. For general n, its form has been 
correctly guessed in Ref. [28] . The oscillatory behavior has also been firstly described 
in Ref. [2] for n — 1, successively its amplitude has been guessed in Ref. [29]. We 
provided an analytical proof of this sequence of numerical guesses and we gave first 
analytical expressions for the oscillating corrections to S n (£) for general n, that was 
too complicated to be guessed. 
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Figure 8. The function F^ at (x) in a semi-infinite chain for different lengths of 
the block £. Finite I corrections are large. 



In Ref. [28], for half-filling and even L, on the basis of numerical results, the 
Renyi entropies have been parametrized as 



Sn 



ln[(L + 1)costtX] + e n - (-l) e b n [L cosirX}- 1/n 



with X = (L/2 - £)/(L + 1). The numerical values b 2 ~ 4.79256 and b 3 ~ 4.19726 
have been reported [28]. Eq. (70) is fully compatible with this expression and predicts 
the value of the amplitude b n for any n. The numerical estimates for b n in Ref. [28] 
coincide with the exact values, showing the high numerical accuracy of Ref. [28]. 
Away from half-filling, in Refs. [28, 29], in order to force a finite-size (or finite-trap) 
scaling, a new parameter has been introduced to take into account the non-scaling 
term sm[k F (2£ + 1)], that is implicit in our definition of hp. The two representations 
are equivalent, but we find Eq. (70) clearer. 



6. Block disconnected from the boundary 

We consider in this section the Renyi entanglement entropies for a block disconnected 
from the boundary. The numerical calculations are straightforward and can be done 
exactly in the same way as before just by considering the correlation matrix C nm 
starting from a spin different from the first. We stress once again that the calculation 
of the spin entanglement can be done simply in terms of fermions, as a peculiarity 
of the chain with open boundary conditions. In the case of two intervals in a PBC 
chain, the entanglement of spins and fermions are known to be different [21, 32] and 
also other boundary conditions are generally expected to make spins and fermions 
inequivalent, so that the calculations should be done following the general method to 
tackle with the Jordan- Wigner string introduced in Ref. [33]. 

For simplicity we will consider only systems in the thermodynamic limit (i.e. 
semi-infinite) and at half filling (hp = tt/2), but the results are very general. In Fig. 
8, we report the function i^ at (a;) obtained dividing Trp^ by the scaling factor in Eq. 
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Figure 9. Scaling of the finite £ corrections for S n (l) vs £ _1 / n for n = 1,2, 3, 5. 
Leading corrections are of the expected form <x and the extrapolation to 

I — > oo clearly gives F n (x) = 1 identically. For large values of n, subleading 
corrections <x l~il n must be included. 



(5), i.e 

c/12(n-l/ri) 



f (*) = — ^^n^-u,^ (71) 



where x is the 4-point ratio in Eq. (6). Global conformal invariance implies that F^x) 
is a function only of x, while in the figure we clearly see different curves for different 
I. As usual [18, 21, 33], these differences are due to finite £ and 1$ effects that are 
severe. In Ref. [25], it has been shown that finite £ corrections are generically of the 
same form (i.e. governed by the same unusual exponent) independently of the number 
of blocks (that in the present case generalizes to its location) . Notice that the data in 
Fig. 8 are not asymptotic also because the various curves do not have the conformal 
symmetry x —> 1 — x [18, 19] (again this is rather common [18, 21, 33] for finite £). 

We can now proceed to the calculation of the asymptotic value of the function 
F n (x) for various n, by using the fact that corrections to the scaling are of the form 
£~ x l n . In Fig. 9, we report for n = 1, 2,3, 5 the function F^ it (x) at fixed x, obtained 
numerically as explained above, as function of £~ 1 / n . For large enough £, the points are 
aligned on straight lines, confirming the correctness of the finite £ scaling (increasing 
n, more corrections of the form £~i/ n must be included to reproduce the numerical 
data, as obvious). This allows to extrapolate to £ — > oo. The result is evident from 
Fig. 9: 

F n (x) = l, (72) 
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identically. At first, this can seem very strange when compared with the complicated 
functions found for two intervals in the XX model with PBC [18, 19, 33]. The simplicity 
of this result is due to the fact that (in the present case) spins and fermions are 
equivalent. For PBC free-fermions, the explicit computation shows F n (x) = 1 [40] 
(as confirmed in some numerical works [41]). This result straightforwardly generalizes 
to OBC free-fermions [42]. The main peculiarity of OBC is not F n (x) — 1, but the 
equivalence of fermions and spins. 

7. Conclusions 

In this manuscript we provided a number of exact results for the asymptotic scaling of 
the Renyi entanglement entropies in open XX spin-chains. Schematically our results 
can be summarized as follows. 

• In the case of a semi-infinite system and a block starting from the boundary, 
we derive rigorously the asymptotic behavior for large block sizes on the basis 
of a recent mathematical theorem for the determinant of Toeplitz plus Hankcl 
matrices. This is given by Eq. (38). 

• To obtain the corrections to the scaling to the asymptotic behavior, we conjecture 
a generalized Fishcr-Hartwig form for these determinants. Eq. (57) gives the 
exact asymptotic behavior of S n (£) at order o(£ _1 ) for any n. Eq. (65) generalizes 
the result for n = oo, i.e. for the largest eigenvalues of the reduced density matrix. 

• By combining these results with conformal field theory arguments, we derive exact 
expressions also in finite chains with open boundary conditions in Eq. (70). 

• In the case of block detached from the boundary, we again use CFT to derive an 
exact expression for the asymptotic S n given by Eq. (5) with F n (x) = 1. 

All these results have been checked against exact numerical computations. 

For infinite chains, in Ref. [7], using a combination of methods based on the 
generalized Fisher-Hartwig conjecture and a recurrence relation connected to the 
Painlcvc VI differential equation, the corrections to S n {t) have been obtained up to 
order 0(£~ 3 ). To the best of our knowledge the recurrence relation obtained in Ref. 
[35] for an infinite system has not been generalized to a semi-infinite one, althought a 
random matrix description of these chains exists [43, 44]. 

It would be also very interesting to consider the exact calculation of asymptotic 
and subleading terms in the Renyi entropies for different boundary conditions and the 
generalization to the XY model on the lines of Refs. [45, 46]. Finally, the accurate 
numerical results of Refs. [28, 29] show a very similar structure also for the case of a 
general confining (trapping) polynomial potential and it is natural to wonder whether 
exact results can be obtained also in those cases. 
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